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For a class of hindered magnetic dipole (Ml) transition processes, such as T(3S I ) — > rji + 7 (the 
discovery channel of the r/f, meson), the emitted photon is rather energetic so that the traditional 
approaches based on multipole expansion may be invalidated. We propose that a "hard-scattering" 
picture, somewhat analogous to the pion electromagnetic form factor at large momentum transfer, 
may be more plausible to describe such types of transition processes. We work out a simple factor- 
ization formula at lowest order in the strong coupling constant, which involves convolution of the 
Schrodinger wave functions of quarkonia with a perturbatively calculable part induced by exchange 
of one semihard gluon between quark and antiquark. This formula, without any freely adjustable 
parameters, is found to agree with the measured rate of T(3S) — > r/t + 7 rather well, and can 
also reasonably account for other recently measured hindered Ml transition rates. The branching 
fractions of T(45) rfc + 7 are also predicted. 
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As a century-old subject, electromagnetic (EM) transi- 
tions have been extensively studied in the fields of atomic, 
nuclear and elementary particle physics. EM transition 
is of considerable interests in heavy quarkonium physics 
from both experimental and theoretical aspects [l[. Ex- 
perimentally, it proves to be a powerful tool to discover 
new quarkonium states that cannot be directly produced 
in e + e~ annihilation into a virtual photon. A very recent 
example is that the long-sought bottomonium ground 
state, the ryt, meson, was finally seen by Babar collab- 
oration in the magnetic dipole (Ml) transition process 
T(35) — > r\hl 0- Theoretically, it provides a useful 
means to probe the internal structure, and the interplay 
between different dynamic scales in quarkonium. 

The standard textbook treatment of EM transitions is 
based on the concept of multipole expansion 1 , by assum- 
ing the emitted photon to be ultrasoft, i.e. k^ 1 ~ mv 2 , 
where m is heavy quark mass, v denotes the typical ve- 
locity of the quark inside a quarkonium. Consequently, 
the long wave length of photon cannot resolve the geo- 
metrical details of quarkonium. Obviously, the multipole 
expansion method is valid provided that kr <C 1, where 
r ~ 1 /mv is the typical radius of a quarkonium. 

One of the great theoretical undertakings is to under- 
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1 The literal meaning of this term is to expand the electromagnetic 
field A M (t , R± jj ) around R, the center-of-mass coordinate of QQ 
pair, in powers of the relative coordinate r, and the expansion 
parameter is essentially kr, k denoting the photon momentum. 
Some authors prefer to dubbing it as long wave-length approxi- 
mation. These two terms are equivalent in this work. 



stand EM transitions in a situation where the multipole 
expansion may break down. It is difficult to find such 
a situation in atomic system, since the typical atomic 
energy spacings are always of order mv 2 ~ ma 2 (where 
a ss 1/137 is the fine structure constant). By contrast, in 
the realm of QCD, the linearly-rising inter-quark strong 
force can host rather highly excited quarkonium states, 
thus one may encounter EM transitions in quarkonium 
with energetic photon. The aim of this work is to offer 
a new perspective to tackle such situation. For definite- 
ness, in this work we will concentrate on the hindered 
magnetic-dipole (Ml) transition {i.e., two quarkonium 
states with the same orbital angular momentum but with 
different spin and principal quantum numbers). Such 
study is of practical importance, because it will help one 
to better understand the process T(3S) — > 77^7, where 
the photon carries a momentum as large as 1 GeV and 
multipole expansion may cease to be a decent method. 

One usually assumes that the Ml transition can pro- 
ceed without gluon exchange between Q and Q. In the 
nonrelativistic limit, the transition rate between two S- 
wave quarkonia is usually described by the well-known 
formula [l[: 
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where eQ is the fractional electric charge of Q, k is the 
photon momentum viewed in the rest frame of n 3 Si 
state, and R n i (r) stands for the radial Schrodinger wave 
function of quarkonium of principal quantum number 
n and orbital angular momentum I. The spherical 
Bessel function jo (jo (2) = ^f 1 ) takes into ac- 
count the so-called finite-size effect (equivalently, resum- 
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ming multipole-exanded magnetic amplitude to all or- 
ders). When k is expected to be ultrasoft, it is then le- 
gitimate to expand this function, and the leading contri- 
bution to the hindered transition vanishes due to orthog- 
onality of wave functions. For several observed hindered 
transition processes, Eq. (JTJ) usually yields predictions a 
few times smaller than the measured values. 

It is widely believed that hindered Ml transitions are 
very sensitive to the relativistic corrections to ([1]) . Unfor- 
tunately, the way of implementing relativistic corrections 
seems to be rather model-dependent. For example, some 
authors proposed that, among the contributions from the 
relativistic corrections to the Ml transition, the hypoth- 
esized scalar-part of the confinement potential may play 
an eminent role, as well as the large anomalous magnetic- 
dipole moment that may be acquired by the bound quark 
due to some nonperturbative mechanism [lj. However, 
both of these suggestions seem not to be based on a firm 
and indisputable footing. As a matter of fact, a vari- 
ety of quite different predictions to the transition rates 
of Y(3S, 2S) — > 77b 7 have been made by different authors 
over the years [3(. When confronting with the recently 
established experimental results, however, most of them 
seem not to be favored. 

Recently, relativistic corrections to the Ml transition 
have been readdressed from the angle of nonrelativistic 
effective field theories (EFT) 4] , which allows one to crit- 
ically examine the validity of some popular, yet maybe 
ad hoc, assumptions in many potential model approaches. 
However, it still remains a great challenge to accommo- 
date the hindered Ml transition in this EFT framework. 
For instance, after including all types of conceivable rel- 
ativistic corrections, the predicted rate for T(2S) —> rjb"f 
seems to be much larger than the measured one. 

Impressive progress has been made in calculating the 
transition rate of J/ip — >■ i] c j directly from lattice QCD 
simulation [5]. However, it is challenging to analyze very 
hindered EM transitions, since excited quarkonium states 
will be difficult to probe by lattice simulation. There is 
also attempt to model the coupled channel effects for 
ip' — > r? c 7, but such framework seems not to be very 
predictive due to existence of several purely phenomeno- 
logical parameters @. 





FIG. 1: Two of four lowest-order diagrams contributing to 
hindered Ml transition in our "hard-scattering" picture. 



In view of shortcomings of the traditional approaches, 
we present a new attempt to analyze very hindered Ml 



transition processes typified by T(3S) — > J]bj- The key 
observation is very simple: in such situation, it is more 
appropriate to count the radiated photon as semihard 



(ki- 



;, often called soft in NRQCD terminology), 



rather than ultrasoft. As a consequence, we should and 
must give up the notion of multipole expansion. We fur- 
ther make a key assumption: the leading contribution to 
such a very hindered transition is described by Fig. [1] 
The underlying rationale is that, in order for the spec- 
tator antiquark to join the final quarkonium state with 
a significant probability, a semihard gluon must be ex- 
changed between Q and Q to exert a kick on it. 

It may be worth digressing into pion EM form factor 
temporarily At large momentum transfer, there exists a 
well-known factorization theorem for this case \7\: 
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where <f>^ implies the nonperturbative light-cone distri- 
bution amplitude of a pion, and T refers to the hard- 
scattering part, which can be computed in perturbation 
theory. The lowest-order contribution to T is also de- 
picted by Fig. Q] It is generally believed that, at large 
Q , this hard-scattering picture is physically more plau- 
sible than the so-called Feynman mechanism (without 
exchange of a hard gluon). 

We plan to derive a factorization formula analogous 
to {2]). In our process, the analogous "hard-scattering" 
part is obtained by integrating out the semihard mode. 
We will assume this part is also perturbatively calcula- 
ble, crucially because rav ^> Aqcd, which seems legiti- 
mate for the T, presumably even for the tp family. Quite 
naturally, we expect that the counterpart of (j) v {x) in 
our nonrelativistic problem will be the Schrodinger wave 
function of quarkonium. 

In passing we highlight the very different role played 
by the semihard mode in this work and in Ref. Q . In the 
latter case when photon is treated as ultrasoft, the semi- 
hard mode can only appear in loop. In contrast to the 
potential mode (p ~ mv 2 , p ~ mi?), it does not make 
contribution when descending from NRQCD onto poten- 
tial NRQCD [4|. According to our scheme, however, the 
semihard mode already makes crucial contribution at tree 
level. It is the very mode that we attempt to integrate 
out perturbatively, in order to fulfill the intended factor- 
ization. 

This said, let us turn to the derivation of the very 
hindered 3 Si — > 1 So radiative transition rate. We will 
perform the calculation in a covariant fashion at the level 
of QCD. Since the hard (p M ~ m) quanta decouple in this 
process, it is also feasible, perhaps more illuminating, to 
directly start from NRQCD. We first note that parity and 
Lorentz invariance constrain the transition amplitude to 
be the form 
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where £[ Tl 35 1 ] and £ 7 represent the polarization vectors 
of the initial quarkonium and the photon, respectively. 
At the rest frame of the initial state, as we will always 
work in, the Lorentz structure becomes e^sg^ • k x e*, 
clearly corresponds to the Ml transition. The scalar co- 
efficient si encodes all the nontrivial dynamics, and we 
will proceed to deduce its explicit form. 

We begin with the parton process Q(p)Q(p) — > 
Q(p')Q(P') + 7(A), as indicated in Fig. [Q We assign the 
momentum carried by each constituent as 
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where q and q' are relative momenta inside each pair, 
which satisfy P ■ q = P' ■ q' = 0. The invariant mass of 
the pairs arc P 2 = 4E 2 and P' 2 = AE 2 , , and the Lorentz 



scalars E q = yr 



q\ Eg, 



q' 2 , which guar- 



antees that each (anti)quark stays on their mass-shell. 
Note in the rest frame of P^ , q^ becomes purely space- 
like. 

The quark propagator in Fig. [TJi) can be expanded: 
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1 



(p' + k) 2 — to 2 
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because k ■ q' ~ m 2 v 2 <C k ■ P' — k ■ P ~ m 2 v. We have 
neglected the small q' component induced by the recoil- 
ing of P' , as well as the Lorentz boost effect on q' , which 
are higher order corrections. The quark propagator in 
Fig. []J>) can be expanded in a similar fashion. Note this 
expansion is also legitimate when k is ultrasoft. 
Fig. QJi) and b) share a common gluon propagator: 



: .(5) 
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Here we retain the ie term explicitly, for the momentum 
integration to be properly evaluated. The two terms in 
the denominator are of comparable size, so ([5]) cannot be 
further expanded. If k is nevertheless counted as ultra- 
soft, the second term can be treated as a perturbation. 
Note our situation is in drastic contrast to the ordinary 
NRQCD calculation for hard exclusive processes. In that 
case, there is always a hard scale > to in the propagators, 
so it is safe to neglect q^'' at the zeroth order of NRQCD 
expansion. 

Having specified the concrete forms of the quark and 
gluon propagators, we then project the quark amplitude 
for Q(p)Q{p) — > Q{p')Q{p') + 7(&) onto the correspond- 
ing color-singlet quarkonium Fock states, with the aid 
of the covariant spin projectors accurate to all orders in 
q(') Q , and include the respective momentum-space wave 
function for each quarkonium state (e.g., see [9(). At the 
lowest order in q and q' , we only need retain the first 



term in (|4|, and neglect all the occurrences of q^ in the 
numerator of the amplitude. It turns out that Fig. [TJi) 
then exactly cancels against Fig. [2)), thus rendering a net 
vanishing result at lowest order in velocity expansion 2 . 

To obtain a nonvanishing prediction, one must proceed 
to the first order in gW in the amplitude. To this level of 
accuracy, it is legitimate to set E q = E q * to, since the 
induced error is of quadratic order in gW. It is a curi- 
ous fact that, if one still keeps only the first term in the 
quark propagator (H]), the O(q^) pieces from the spin 
projectors and the ^ term from the quark propagator 
then make a nonzero contribution in an individual dia- 
gram, but their contributions again cancel upon summing 
Fig. QJi) and Fig. [TJ>). Therefore, the leading surviving 
0(q^) contribution can only be obtained by retaining 
the second term in the expanded quark propagator ([5]), 
while neglecting gW terms altogether in elsewhere of the 
amplitude. After some efforts, we can read off the re- 
duced amplitude: 
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where Cp = |, and the prefactor 2 indicates that two 
undrawn diagrams make equal contributions as Fig.Q}i) + 
b), owing to charge conjugation symmetry. </>„</) o signifies 
the momentum-space Schrodinger wave function, and the 
"hard-scattering" kernel is 



kq 



it 



(7) 



Eq. (|6|) is the desired factorization formula in momentum 
space. 

It would be more convenient to work with the famil- 
iar spatial wave functions. Thanks to the fact that the 
"hard-scattering" part depends only on the difference be- 
tween the relative momenta of two quarkonia, q' — q, and 
not on q or q' separately, upon Fourier transformations, 
one can arrive at a compact expression in the position 
space via contour integral: 



_ Aee Q C F a s 
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(8a) 



drr 2 R* n , (r),?(r)R n0 (r), (8b) 



where R(r) appearing in the overlap integral S nn i is the 
radial wave function. We have used the relation k ■ P = 
k M n and M n is the mass of the initial-state quarkonium. 
The dimensionless kernel 2?(t) is obtained by Fourier 



2 This is somewhat analogous to the hard exclusive process rji, — > 
J/tpJ/tp, where the amplitude also vanishes at the lowest order 



charm quark relative velocity Q. 
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transforming T(q) and integrating over solid angle: 

. ei kr [ fkr\ 2 / kr\ ( kr\ 

W = M7r[ JO {Y)-Vr Jl {Y) +lJl {Y) ' 

(9) 

where ji(x) = ^£ — £2*2., jt may be worth remind- 
ing that the above combination of spherical Bessel func- 
tions in the bracket resembles the conventional electric- 
dipole (El) transition formula with finite-size effect in- 
corporated. Notice S?(r) develops an imaginary part, 
since the exchanged semihard gluon can become on-shell 
when q q' = k. However, we would like to stress that, 
the characteristic virtuality of the exchanged gluon in fl5J 
should be of order m 2 v 2 3> Aq CD , thus the emergence of 
imaginary part in the "hard-scattering" kernel should be 
viewed as an artifact due to ignoring the recoiling effect of 
P'. If the effect of the imaginary part is insignificant with 
respect to that of the real part, we may feel such a igno- 
rance is tolerable, otherwise it will indicate a theoretical 
disaster. As we will see in later phenomenological anal- 
ysis, the contamination of the imaginary part is indeed 
always negligible for a class of hindered Ml transitions 
in bottomonium and charmonium systems. 

It is interesting to examine the asymptotic form of 
3T(r) as kr <C 1. Using ji(x) ~ ( 2 ;+i)!! a ^ sma U x > 
one finds 2?(r) —> + i-^ as kr —> 0. It turns out 
that, the real part might be identified with, up to a con- 
stant, the O(o. s ) matching coefficient y^ T '( rxrx - B )]/™ 1 m 
Sec. IIIC of [4J (note there arises some subtle issue regard- 
ing gauge invariance) . Since the imaginary part becomes 
r-independent in the long wavelength limit, as expected, 
it does not contribute to the hindered Ml transition. 

Finally we can express the transition width as 

r[n 3 5 1 ^n /1 5 + 7 ] = 4-K| 2 

LZTT 

= f oe^CjaJK*,! 2 , (10) 

where we have averaged upon spin of the initial 3 Si state 
and sum over two transverse polarizations of the photon. 

Eq. (|10J is the key formula of this work, which 
looks quite simple. In evaluating the overlap integral 
S nn i i the input wave functions are obtained by solving 
Schrodinger e quat ion with the widely-used Cornell po- 
tential model [10| and Buchmuller-Tye (BT) potential 
model 11]. Parameters in both potential models are 
tuned such that the bb and cc spectroscopy below open 
flavor threshold are successfully reproduced. The only 
freely adjustable parameter seems to be the strong cou- 
pling constant, a s (fi). However, the choice of the renor- 
malization scale [x is by no means arbitrary. On physical 
ground, it should be fixed around the typical value of 
the quark 3-momentum in quarkonium, which is about 
1.2 GeV for bb system, and 0.9 GeV for cc system [l2T |. 
Therefore, with a s fixed, our formalism becomes rather 
predictive, and, readily falsifiable. 



In TableUwe have tabulated various predictions to hin- 
dered Ml transitions of n 3 Si —> n n So- We also present 
the numerical results for $ nn i , and reassuringly, the con- 
tribution from Im^(r) is indeed insignificant. As one can 
tell, the agreement between our predictions, especially 
from the Cornell potential model, and the measurement 
for the transition rate of Y(3S) — > rjbj, is strikingly suc- 
cessful. Curiously, for other hindered Ml transitions, 
where the photon is not that energetic so that the mul- 
tipole expansion method may still apply, our formalism 
again appears to make a decent account of the measured 
transition rates, agrees typically within 2 — 3 a. It seems 
fair to conclude that our simple factorization formula has 
passed quite nontrivial tests. Given the fact that there 
is almost no free parameters in (|10j) . we feel encouraged 
that our formalism has captured at least some correct 
and relevant ingredients. We hope future measurements 
of T(4S) — > 77f,7 can further test our mechanism. 

It might be tempting to seek simplified expression for 
the overlap integral £" n n', by exploiting some hierarchy 
between different bb energy levels. Higher radial excita- 
tion, say, T(3S), is known to have considerably larger 
radius than rji,. An intuitive guess is that (031 may not 
be necessarily sensitive to the full profile of Rso(r), in- 
stead may only sensitive to its value at small distance 
(about the radius of rib). If this were true, one could pull 
R^o( r ) outside of the integral, and approximate it by its 
value at the origin. The transition rate predicted this way 
turns out to be about two orders of magnitude greater 
than the measured one! If we play the same game for 
Rl (r), the result would be about ten times larger than 
the data. The failure of these approximations may be 
understood from the empirical fact that, in the Cornell 
or BT potential models, the average momentum of quark 
in different bb energy levels is more or less equal. As a 
result, there seems no ground to neglect q or q' in the 
"hard-scattering" kernel in ©. 

For the Ml transition from n So to n' 3 Si, one needs 
multiply (HHJ) by a statistical factor of 3. Various par- 
tial widths for rjb(nS) — > T7 are about 10 eV, and that 
for rj c (2S) — > J/tp'f is about 1 keV. These bottomonium 
transitions may be accessible in high-energy hadron col- 
lider experiments such as CERN Large Hadron Collider 
(LHC), and BESIII program may provide a chance to 
look for this charmonium hindered transition. 

As in any factorization framework, we expect that 
the factorization formula ^ is perturbatively improv- 
able. It will be a major progress to calculate the next-to- 
leading order correction to the "hard-scattering" kernel. 
To achieve this, it might prove easier to reformulate our 
derivation in the context of NRQCD. It would also be 
interesting to implement relativistic corrections to ([6]). 

Obviously our strategy needs not to be confined to hin- 
dered Ml transitions only. It should be applicable when- 
ever the radiated photon cannot be viewed as ultrasoft 
and multipole expansion breaks down. It will be interest- 
ing to work out the corresponding factorization formula 
for El transitions such as Xbj(2P) —> T7. It would be 
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TABLE I: Measured and predicted branching fractions of various hindered Ml transition processes n 3 Si — > n' So + 7 for 
bottomonium and charmonium. The photon momentum k is determined by physical kinematics. The total widths of T(nS) 
and tp(2S) states, as well as all the quarkonium masses, are taken from PDG08 compilation [13J, except rj b mass is taken to 
be 9389 MeV @], and rj b (2S) mass taken as 9997 MeVjjj]. For T(2S)-fyn b , we use the preliminary Babar result 0; for 
ip{2S)^ r )Tj c , we quote the latest Cleo measurement [15|. instead of the world average value given in [l3|]. We have taken 
a a (fJ,) = 0.43 and 0.59 for fj, = 1.2 and 0.9 GeV, respectively. 
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B (Exp.) 


a s 




(xlO" 2 ) 


B ( Our predictions) 


modes 


(MeV) 






Cornell 


BT 


Cornell 


BT 


T(2S , )-> 7 r ?6 


614 


(4.2 ± 1.4) x 10" 4 


0.43 


3.7e i2 °" 


3.2e i2 - 7 " 


1.4 x 10~ 4 


1.1 x 10 -4 


T (35)^,6 


921 


(4.8 ± 1.3) x 10~ 4 


0.43 


27e i2.e° 


2.3e a5 ° 


3.7 x 10" 4 


2.8 x 10" 4 




1123 




0.43 


2.2e !2 ' 8 ° 


1.9e l3 - 7 ° 


4.3 x 10~ 7 


3.2 x 10" 7 


T(4S)^77 7b (2S) 


566 




0.43 


1.7e !2 - 2 ° 


1.6e i2 ' 7 ° 


3.2 x 10" 8 


2.7 x 10" 8 


^(25 , )->77 7c 


638 


(4.3 ±0.6) x 10" 3 


0.59 


6.4e l9 - 7U 


5.7e' 12 ' 9U 


2.7 x 10~ 3 


2.1 x 10~ 3 



also interesting to generalize this "hard-scattering" for- 
malism to explore the hadronic transition processes such 
as T(3S*, 4S*) T + nn. 
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